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Abstract. The multipole analysis investigates the 
arrival directions of registered neutrino events in 
AMANDA-II by a spherical harmonics expansion. 
The expansion of the expected atmospheric neutrino 
distribution returns a characteristic set of expansion 
coefficients. This characteristic spectrum of expan- 
sion coefficients can be compared with the expansion 
coefficients of the experimental data. As atmospheric 
neutrinos are the dominant background of the search 
for extraterrestrial neutrinos, the agreement of ex- 
perimental data and the atmospheric prediction can 
give evidence for physical neutrino sources or sys- 
tematic uncertainties of the detector. Astrophysical 
neutrino signals were simulated and it was shown 
that they influence the expansion coefficients in a 
characteristic way. Those simulations are used to 
analyze deviations between experimental data and 
Monte Carlo simulations with regard to potential 
physical reasons. The analysis method was appUed on 
the AMANDA-II neutrino sample measured between 
2000 and 2006 and results are presented. 

Keywords: Neutrino astrophysics, Anisotropy, 
AMANDA-II 

I. Introduction 

The AMANDA-II neutrino detector located at South 
Pole was constructed to search for astrophysical neutri- 
nos. These neutrinos could originate from many different 
Galactic and extragalactic candidate source types such 
as Active Galactic Nuclei (AGN), supernova remnants 
and microquasars. The detection of neutrinos is based on 
the observation of Cherenkov light emitted by secondary 
muons produced in charged current neutrino interactions. 
This Ught is observed by photomultipliers deployed in 
the Antarctic ice. Their signals are used to reconstruct 
the direction and the energy of the primary neutrino. 

AMANDA-II took data between 2000 and 2006. The 
background of atmospheric muons is reduced by select- 
ing only upward-going tracks in the detector, as only 
neutrinos are able to enter the detector from below. This 
restricts the field of view to the northern hemisphere. 

The data is filtered and processed to reject 
misreconstructed downward-going muon tracks |1 1. The 
final data sample contains 6144 neutrino induced events 
between a declination of 0° and +90° with a purity of 
> 95% away from the horizon. 



II. Analysis Principle 

The idea of this analysis is to search for deviations of 
the measured AMANDA-II neutrino sky map from the 
expected event distribution for atmospheric neutrinos, 
which constitute the main part of the data sample fj). 
A method to study such anisotropics is a multipole 
analysis, which was also used to quantify the Cosmic 
Microwave Background fluctuations. The analysis is 
based on the decomposition of an event distribution 
fi^.^t^) = E»=T"'' '^(cose*,; - cos6')(5(0, - 0) into 
spherical harmonics Y;™(0, 0), where 6 and are the 
zenith and azimuth of the spherical analysis coordinate 
system. The expansion coefficients are 

aj" = / d(t) I dcose f{e,4>)Yi"'*{e,(f>). (i) 

Jo J-1 

They provide information about the angular structure of 
the event distribution f{9,(f)). The index I corresponds 
to the scale of the angular structure 5 k, while 
m gives the orientation on the sphere. The expansion 
coefficients with m = depend only on the structure in 
the zenith direction of the analysis coordinate system. 
Averaging over the orientation dependent a™ yields the 
multipole moments 

+' 

Ci = V la^l^. (2) 

m— — l 

They form an angular power spectrum characteristic for 
different input neutrino event distributions. 

The initial point of this analysis is the angular power 
spectrum of only atmospheric neutrino events. There- 
fore, neutrino sky maps containing 6144 atmospheric 
neutrino events according to the Bartol atmospheric 
neutrino flux model |3| are simulated and numerically 
decomposed with the software package GLESP ID. 
Statistical fluctuations are considered by averaging over 
1000 random sky maps, resulting in a mean (C;) and a 
statistical spread of each multipole moment. 

The same procedure is applied to simulated sky maps 
containing atmospheric and different amounts of signal 
neutrinos with a total event number of likewise 6144 
events. The influence of the signal neutrinos on the 
angular power spectrum is studied in terms of the pulls 

J i^l) ~ (C'i,atms) 

di = . {i) 
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1. (a): Pull plot for the multipole moments C; of the isotropic point source model. Sources are simulated with a mean source strength 



Fig. 

fj, = 5 and an EiJ'^ energy spectrum. The number of sources Afgourccs on the full sphere is varied. The corresponding number of signal 
neutrinos on the northern hemisphere is given in brackets. The errors bars are hidden by the marker symbols, (b): Pull plot for the expansion 
coefficients of the cosmic ray interaction model with the Galactic plane in Galactic coordinates. The fraction of neutrinos in the sky map 
originating from the Galactic plane is varied. The corresponding number of signal neutrinos is given in brackets. The errors bars are hidden by 
the marker symbols. 



III. Signal simulation 

The different models for candidate neutrino sources 
investigated in this analysis are: 

1) Isotropically distributed point sources 

2) A diffuse flux from FR-II galaxies and blazars |5 | 

3) AGN registered in the Veron-Cetty and Veron 
(VCV) catalog IS 

4) Galactic point sources such as supernova remnants 
or microquasars 

5) Cosmic rays interacting in the Galactic plane. 

All simulated pointlike neutrino sources are character- 
ized by a Poissonian distributed source strength with 
mean /i and an energy spectrum i?^^. The relative angu- 
lar detector acceptance depends on the neutrino energy 
and therefore on the spectral index of the simulated neu- 
trino source. Signal neutrinos are simulated according to 
this acceptance considering systematic fluctuations. The 
total number of signal neutrinos in a sky map of the 
northern hemisphere with A'sourcos simulated sources on 
the full-sky is therefore given by 0.5 • /z • iVgourcos- 
Additionally the angular resolution is taken into account. 
It dominates over the uncertainty between the neutrino 
and muon direction. 

The spectral index of pointlike sources is varied 
between 1.5 < 7 < 2.3. As the spectral index of 
atmospheric neutrinos is close to 3.7, signal and back- 
round neutrinos underlie different angular detector ac- 
ceptances. Thus, additionally to the clustering of events 
around the source directions also the shape of the total 
angular event distribution is used to identify a signature 
of signal neutrinos in the angular power spectrum |7|. 

Neutrinos from our Galaxy disturb the atmospheric 
event distribution by their bunching within the Galactic 
plane modeled by a Gaussian band along the Galactic 
equator. Neutrinos produced in cosmic ray interactions 



with the interstellar medium of our Galaxy are assumed 
to follow the iJ^^^ primary energy spectrum. 

A further topic (model 6) that can be studied with a 
multipole analysis are neutrino oscillations. The survival 
probability of atmospheric muon neutrinos depends 
on the neutrino energy and the traveling length of the 
neutrino as well as the mixing angle 623 and the squared 
mass difference Atojs- The traveling length can be 
expressed by the Earth's radius and the zenith angle of 
the neutrino direction [7|. Thus, the neutrino oscillations 
disturb the angular event distribution of atmospheric 
neutrinos. With the assumption of sin^ (26'23) ~ 1 the 
squared mass difference remains for investigation. Due 
to the relatively high energy threshold of 50 GeV the 
effect is small. 

IV. Evaluation of the power spectra 

The deviations from a pure atmospheric angular power 
spectrum caused by signal neutrinos are studied by the 
pulls. These pulls are exemplarily shown in Fig. [Talfor 
the model of isotropic point sources. The behaviour 
of the pulls is characteristic for each signal model. 
Different multipole moments carry different sensitivity 
to the neutrino signal. The absolute value of the pull 
increases linearly with the amount of signal neutrinos in 
the sky maps. Each pull has a predefined sign. 

The deviation of a particular sky map with multipole 
moments Ci from the pure atmospheric expectation 
(G.atms) is quantified by a significance indicator 
defined as 

7-i2 1 f ~ {Cl.a.tms)\ 

D =-, 2^ sgn; -wi-l ■ , (4) 

'max \ "'Cl.atma / 

where Zmax determines the considered multipole mo- 
ments. The term in brackets is the pull between the 
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particular sky map and the mean of the atmospheric 
expectation as defined in Eq. [3] The factors 



Wl = 



(Ci) — (C;^atms) 



(5) 



are defined to weight the pulls according to their ex- 
pected sensitivity to the signal. For each neutrino signal 
model one dedicated set of weights wi is determined. 
Due to the linear increase of the pulls with the signal 
strength the strength chosen to calculate wi is arbitrary. 

The weight factors wi carry the expected sign of the 
pulls, sgn; is the sign of the measured pull. Thus, the 
calculated for the particular sky map is increased if 
the observed deviation has the direction expected for the 
signal model and reduced otherwise. 

Due to the weighting of the pulls, the sensitivity 
becomes stable for high Zmax- A choice of Zmax = 100 
is sufficient to provide best sensitivity to all investigated 
signal models. 

The of a sky map is interpreted physically by the 
use of confidence belts. Therefore, 1000 sky maps for 
every signal strength within a certain range are simu- 
lated and the D^-value for each sky map is calculated 
separately to obtain the distributions. The calculation 
of the average upper limit at 90% confidence level 
assuming zero-signal is used to estimate the sensitivity 
of the analysis to different astrophysical models apriori. 

As the multipole analysis is applied to a wide range 
of astrophysical topics, the trial factor of the analysis 
becomes important. The trial factor raises with each new 
set of weights used to evaluate the experimental data. 
For this reason, models with almost similar weights are 
combined to a common set of weights and only six sets 
are remaining. 

If the signal signatures show up only in the zenith 
direction of the analysis coordinate system the expansion 
coefficients are more sensitive than the multipole 
moments C/. The reason is, that the expansion coeffi- 
cients with m = are independent from the azimuth 
(j) and contain the pure information about the zenith 
direction 9. A signal only depending on 6 causes only 
statistical fluctuations but no physical information in 
the other expansion coefficients. Therefore, the signal 
has only power in the a°. The analysis method stays 
exactly the same in these cases, except that all C/ are 
replaced by the a°. This is related to the models of 
neutrinos from the Galactic plane and from sources of 
the VCV catalog, which show north-south-symmetries 
of the neutrino signals in Galactic and supergalactic 
coordinates, respectively. Unlike the multipole moments 
Ci, the a° do not average over different orientations. 
Therefore, the analysis of the af strongly depends on the 
used coordinate system. An example for pulls of a° for 
the model of a diffuse neutrino flux from the Galactic 
plane is shown in fig. [lb] The characteristic periodic 
behavior of the pulls is explained by the symmetry 
properties of the spherical harmonics. 
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Fig. 2. Pull plot for the experimental multipole moments Ci . Expected 
pulls for typical model parameters of isotropic point sources ai'e shown 
for comparison. The error bars symbolize the statistical fluctuation 
expected for an atmospheric neutrino sky map. 
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Fig. 3. Pull plot for the experimental expansion coefficients in 
Galactic coordinates. Expected pulls for typical parameters of cosmic 
ray interactions with the Galactic plane are shown for comparison. 
The error bars symbolize the statistical fluctuation expected for an 
atmospheric neutrino sky map. 



V. Experimental results 

The experimental data is analyzed in two steps. First, 
the experimental data is tested for its compatibility with 
the pure atmospheric neutrino hypothesis. Secondly, the 
experimental pulls are compared with the expectations 
for the different investigated neutrino models. 

The pulls of the experimental data are shown for the 
multipole moments Ci in Fig. |2] and for the expansion 
coefficients a" in Galactic coordinates in Fig. |3] To 
compare the measured data with the expected event dis- 
tribution, a is calculated for the multipole moments 
Ci and the expansion coefficients a" for transformations 
into equatorial. Galactic and supergalactic coordinates 
separately. As no signal model is tested sgn^ = = 1 
is assumed. A comparison with the corresponding 
distributions results in the p-values giving the proba- 
bility to obtain a which is at least as extreme as 
the measured one assuming that the pure atmospheric 
neutrino hypothesis is true (Table IJl. 
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The statistical consistency of C; and in equa- 
torial coordinates with the atmospheric expectation is 
marginal. Rotating to inclined coordinate systems, e.g. 
Galactic and supergalactic, the consistency improves. 
The deviation from the pure atmospheric expectation is 
not compatible with any of the signal models (see Fig. 
|2] [3] for examples). The discrepancy may be attributed 
to uncertainties in the theoretical description of the 
atmospheric neutrino distribution, or to a contribution 
of unsimulated background of down-going muons mis- 
reconstructed as up-going, or to the modeling of prop- 
erties of the AMANDA detector. 

TABLE I 

P-VALUES FOR THE COMPATIBILITY OF EXPERIMENTAL DATA AND 
PURE ATMOSPHERIC NEUTRINO HYPOTHESIS. 



Observable 


p- value 


Ci 

a9. Equatorial 
a?. Galactic 
ay, Supergalactic 


0.02 
0.02 
0.15 
0.70 



The signal models are tested by calculating the De- 
values of the experimental data using the corresponding 
sign and weight factors. As the observed deviations do 
not fit any of the investigated signal models the physical 
model parameters are constrained. Due to the observed 
systematic effects affecting mainly the multipole mo- 
ments Ci and the equatorial expansion coefficients a'} 
no limits on the models analyzed in the corresponding 
coordinate systems (models 1, 2 and 6) are derived. The 
other models are less affected. The limits given below 
do not include these systematic effects. 

A limit on the source strength assuming the VCV 
source distribution (model 3) is calculated for those 
sources closer than 100 Mpc to the Earth. In this model 
all sources are expected to have the same strength and 
energy spectrum. For a typical spectral index of 7 = 2 
the average source flux is limited by the experimental 
data to a differential source flux of d^/dE ■ E"^ < 1.6 • 
10~^° GeVcm"^ s^^ sr^^ in the energy range between 
1.6 TeV and 1.7 PeV. 

For the random Galactic sources (model 4), the 
number of sources is constrained assuming the same 
source strength and energy spectrum for all sources 
as well. For a spectral index of 7 = 2, the limit 
on the number of sources is set by AMANDA to 
^sources < 39 assuming a source strength of d^/dE ■ 
E^ < 10-8 GeVcm-2 g-i si-^ or A^sourcos < 4300 for 
sources with d'P/dE ■ E'^ < 10'^° GeVcm-^ g-i sr^i. 
For source fluxes in between the limit can be ap- 
proximated by assuming linearity between A'sourccs and 
log{d^/dE ■ E^). 

The differential flux limit obtained from the 
experimental data on the diffuse neutrino flux from 
cosmic ray interactions in the Galactic plane (model 5) 
is d^/dE ■ E^-^ < 3.2 ■ 10"'' GcV^'^ cm^^ g-i sr^i. 
This flux limit is shown in Fig. |4] together with the 
results of two other AMANDA analyses and two 
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Fig. 4. Limit of the 7yr multipole analysis on the diffuse neutrino flux 
from cosmic ray interactions in the Galactic plane in dependence of 
the valid energy range. The limit is compared with two other analyses 
|2|, |8 | and two theoretical predictions 0, ITOl . 

theoretical flux predictions. The seven year multipole 
analysis provides currently the best limit. However, it 
is still not in reach of the theoretical predictions. 

VI. Conclusion 

It is shown that the multipole analysis is sensitive to 
a wide range of physical topics. Its area of application 
is in particular the field of many weak sources in 
transition to diffuse fluxes. With the statistics of seven 
years of AMANDA data and improvements of the 
analysis technique the method is now restricted by 
systematic uncertainties in the atmospheric neutrino 
zenith distribution of the order of a few percent. 
Transforming to coordinate systems less affected by 
the equatorial zenith angle such as the Galactic and 
supergalactic system physical conclusions are still 
possible. A compatibility of the measurement with 
the background expectation of atmospheric neutrinos 
is observed. Current efforts to better understand the 
observed systematics would allow an application of the 
multipole analysis on future high statistic IceCube data. 
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